TVibology of the lubricant quantized-sliding state 
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In the framework of Langevin dynamics, we demonstrate clear evidence of the peculiar quantized 
sliding state, previously found in a simple ID boundary lubricated model [Phys. Rev. Lett. 97, 
056101 (2006)], for a substantially less idealized 2D description of a confined multi-layer solid lubri- 
cant under shear. This dynamical state, marked by a nontrivial "quantized" ratio of the averaged 
lubricant center-of-mass velocity to the externally imposed sliding speed, is recovered, and shown to 
be robust against the effects of thermal fluctuations, quenched disorder in the confining substrates, 
and over a wide range of loading forces. The lubricant softness, setting the width of the propagating 
solitonic structures, is found to play a major role in promoting in-registry commensurate regions 
beneficial to this quantized sliding. By evaluating the force instantaneously exerted on the top plate, 
we find that this quantized sliding represents a dynamical "pinned" state, characterized by signif- 
icantly low values of the kinetic friction. While the quantized sliding occurs due to solitons being 
driven gently, the transition to ordinary unpinned sliding regimes can involve lubricant melting due 
to large shear-induced Joule heating, for example at large speed. 

PACS numbers: 68.35.Af, 05.45.Yv, 62.20.Qp, 81.40.Pq, 46.55.+d 



I. INTRODUCTION 

When a confined lubricant is characterized by a large 
surface-to-volume ratio, the surface features of the con- 
fining walls may strongly affect its tribological proper- 
ties. If the width of the lubricant film is reduced to a 
few atomic layers, as in the boundary-sliding regime, the 
atoms in the film tend generally to be ordered into layers 
parallel to the bounding substrates^ 2 - Both numerical 
simulations and experiments generally conclude that, in 
such a strongly confined geometry, the film behaves like 
a solid, even at temperatures significantly higher than its 
bulk melting temperature^ In real operative conditions, 
moreover, the case of "dry" friction is quite exceptional. 
A physical contact between two solids is generally me- 
diated by so-called "third bodies" , which act like a lu- 
bricant filmi^ For crystals sliding on crystals, one may 
consider the moving contact as characterized schemati- 
cally by three inherent length scales: the periods of the 
bottom and top substrates, and the period of the embed- 
ded lubricant structure. 

The problem of boundary lubricated friction is fasci- 
nating both from the fundamental point of view and for 
applications. Interesting dynamical behaviors, with pos- 
sible tribologically important implications of an irregu- 
lar distribution of the lubricant velocity in between the 
sliding surfaces, have recently been observed in numeri- 
cal simulations, depending on the "degree" of geometri- 
cal incommensurability defining the moving interfaced— 
The prominent nontrivial feature of those simulation 
is an asymmetry in the relative sliding velocity of the 
intermediate lubricating sheet relative to the two sub- 



strates. Strikingly, this velocity asymmetry takes an ex- 
actly quantized value which is uniquely determined by 
the incommensurability ratios and is insensitive to all 
other physical parameters of the model. The occurrence 
of this surprising and robust regime of motion, giving 
rise to perfectly flat plateaus in the ratio of the time- 
averaged lubricant center-of-mass (CM) velocity to the 
externally imposed relative speed, was ascribed to the 
intrinsic topological nature of this quantized dynamics. 
The phenomenon, investigated in detail in a rather ideal- 
ized ID geometry^— was explained by the corrugation 
of a sliding confining wall rigidly dragging the topologi- 
cal solitons (kinks or antikinks) that the embedded lubri- 
cant structure forms with the other substrate. Evidence 
of the existence of this peculiar regime of motion was 
then confirmed shortly after for a substantially less ide- 
alized 2D model of boundary lubrication, where atoms 
were allowed to move perpendicularly to the sliding di- 
rection and interacted via LJ potentials 15 . The solitonic 
transverse corrugations of the lubricant, propagating in- 
side the film from bottom to top, were observed to favor 
the kink (or antikink) tendency to pin to the top-layer 
spatial periodicity, thus strengthening the quantization 
mechanism. 

In this work, we investigate the tribological signa- 
ture of the quantized- velocity state for such a 2D multi- 
layer solid lubricant under shear. Its robustness against 
thermal effects (implemented by means of a standard 
Langevin approach), quenched disorder in the confining 
substrates, and over a wide range of loading forces is ana- 
lyzed in detail. The lubricant softness, setting the width 
of the propagating solitonic structures, together with the 
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FIG. 1: (color online). A sketch of the model with the rigid 
top (solid circles) and bottom (open) crystalline sliders (of 
lattice spacing a t and ab respectively), the former moving 
at externally imposed a;- velocity « e xt- One or more solid lu- 
bricant layers (shadowed) of rest equilibrium spacing ao are 
confined in between. 



soliton coverage number, is considered in promoting in- 
registry commensurate regions beneficial to the quantized 
sliding. Hysteretic termination and kinetic friction are 
properly evaluated to highlight tribological differences 
with respect to ordinary unpinned sliding regimes. 

The paper is organized as follows. In Sec. UH we 
summarize the main features of the 2D-implemented tri- 
bological confined model, together with the numerical 
procedure adopted to solve the system of coupled equa- 
tions describing the dissipative particle dynamics and top 
plate motion at finite temperatures. Section IIIII is de- 
voted to prove the robustness of the quantized dynamical 
regime against many physical parameters of the model. 
In Sec. IIVI the relative lubricant-lubricant and lubricant- 
substrate interactions are varied to study the role of the 
confined film softness. By evaluating the force instan- 
taneously exerted on the advancing plate in Sec. [Vj we 
characterize the tribological features of the quantized dy- 
namics and show that it corresponds to significantly low 
kinetic friction force. The transition towards ordinary 
sliding regimes can involve lubricant melting due to large 
shear-induced Joule heating, leading to a higher-friction 
sliding state. As detailed in Sec. IVI1 within the frame- 
work of our rigid, constant-speed driving, we observe an 
intermittent stick-slip dynamics, associated to tiny am- 
plitude fluctuations, in the kinetic friction force experi- 
enced by the top substrate while in the quantized- velocity 
state. Discussions and conclusions are given in Sec. IVIH 



II. THE 2D CONFINED LUBRICATED MODEL 

Our system is composed of two rigid walls made of 
equally-spaced atoms and N\ identical lubricant atoms 
confined in between, organized in A^ ayer layers, see Fig.Q] 
where N\ aycl = 5. The lubricant atoms move under the 
action of pairwise Lennard- Jones (LJ) potentials 
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describing the reciprocal interactions. The cutoff radius 
is set at r = r c = 2.49 a a , where $ a (r c ) ~ -8.4 • 10~ 3 e a . 

For the two substrates and the lubricant we assume 
three different kinds of atoms, and characterize their mu- 
tual interactions ($bb ®n an d $ti refer to potential en- 
ergies for the bottom- lubricant, lubricant-lubricant, and 
top-lubricant interactions, respectively) with the follow- 
ing LJ radii a a 



cti = at 



f bi = a b , and (Tii = Q-0 ■ 



(2) 



which, for simplicity, are set to coincide with the fixed 
spacings a t and ab between neighboring substrate atoms, 
and the average x-separation ao of two neighboring lubri- 
cant atoms, respectively. This restriction is only a matter 
of convenience, and is not essential to the physics we are 
describing. These three different periodicities a t , ao and 
ab define two independent ratios: 
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the latter of which we take closer to unity, 
max(Ab,A b _1 ) < A t , so that the lubricant is less 
off-register to the bottom substrate than to the top. 

For simplicity, we fix the same LJ interaction energy 
e t i = en = Cbi = e for all pairwise coupling terms and the 
same mass m of all particles. We take e, ao, and m as 
energy, length, and mass units. This choice defines a set 
of "natural" model units for all physical quantities: for 
instance velocities are measured in units of e 1 / 2 to -1 / 2 . 
In the following, all mechanical quantities are expressed 
implicitly in the respective model units. 

The interaction with the other lubricant and sliders' 
particles produces a total force acting on the j-th lubri- 
cant particle, of coordinate fj, given by 



F 3 = 
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where fa and fbj are the positions of the N% top and 
Ab bottom atoms. By convention, we select the frame 
of reference where the bottom layer is static, with the 
atoms located at 



rbix(t)=ia h , r hiz (t)=0. 



(5) 



The top slider is rigidly driven at a constant horizontal 
velocity i> cx t, and can also move vertically (its inertia 
equals the total mass N t m of its atoms) under the action 
of the external loading force F applied to each particle 
in that layer and that due to the interaction with the 
particles in the lubricant layer: 



rtix(t) = sot 



t. 



mz(t) =r tz (t), (6) 
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where the equation governing r t z is 



N t m f t ■ 
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-N t F. (7) 



To remove the Joule heat, and control the lubricant 
temperature in this driven system, rather than a Nose- 
Hoover thermosta t 16 i 17 as in Ref. [[a, we use a standard 
implementation of the Langevin dynamics, 1 - with the ad- 
dition of a damping term plus a Gaussian random force 
fj(t) to the Newton equations for the lubricant particles. 
The damping force includes the contributions of the en- 
ergy dissipation into both individual substrates: 



/damp i = -If] fi - T) (fi - f t ) . 



(8) 



Taking into account this double contribution to the rj- 
dissipation, the Gaussian null-average random forces sat- 
isfy the relation 



(f jx (t)f jx (t'))=4 V k B TS(t-t') 



(9) 



and similarly for the z components, such that in a static 
configuration (i> e xt — 0) the Langevin thermostat leads 
to a steady state characterized by standard Boltzmann 
equilibrium 



(E k ) =2N^k B T, 



(10) 



with the total kinetic energy of the lubricant E k = 

2 TO ^i • This algorithm represents a simple but nu- 
merically stable and effective phenomcnological approach 
to describe energy dissipation into the substrates occur- 
ring through the excitation of phonons and (in the case 
of metals) of electron-hole pairs. As we shall see below, 
the Langevin thermostat controls and compensates rea- 
sonably well the expected Joule heating of the substrate, 
even in the sliding regime (v cxt ^ 0). 

We adopt a fairly small value 77 = 0.1 leading to un- 
derdamped lubricant dynamics. To guarantee a detailed 
force balance (Newton's third law) , we add the following 
force term 



Ni 



ft) =r/Ni (v cm - f t ) 



(11) 



to the equations for the motion of the top layer. While 
the z-component of this additional term has a real influ- 
ence on the motion governed by Eq. ([7]), of course its x- 
component only affects the external force F k required to 
maintain the velocity component f t x constant and equal 

tO V cxt . 

We integrate the ensuing equations of motion within a 
x-periodic box of size L = N\ao, by means of a standard 
fourth-order Runge-Kutta method, 19 We usually start off 
the dynamics of a single lubricant layer from equally- 
spaced lubricant particles at height i\ z = ab and with the 



top layer at height r tz = ab + a t , but we consider also dif- 
ferent initial conditions. For several layers, we initialize 
the system with lubricant particles at perfect triangular 
lattice sites, and the top slider correspondingly raised. 
After an initial transient, sometimes extending for sev- 
eral hundred time units, the sliding system reaches its 
dynamical steady state. In the numerical simulations, 
an adiabatic variation of the external driving velocity 
(a controllable parameter in tribological experiments) is 
considered and realized by changing w cxt in small steps 
5v = 0.1, letting the system evolve at each step for a time 
long enough for all transient stresses to relax, in practice 
1500 time units. In the calculations presented below, we 
compute accurate time-averages of the physical quanti- 
ties of interest by averaging over a simulation time of 
3000 time units or more after the transient is over. At 
higher temperature, fluctuations of all physical quantities 
around their mean values increase, sometimes requiring 
even longer simulation times to obtain well-converged av- 
erages. To estimate error bars of all average quantities, 
we split the whole trajectory into 30 pieces of equal du- 
ration, and then evaluate the standard deviation of the 
averages carried out over each individual interval. 



III. THE ROBUSTNESS OF THE PLATEAU 
DYNAMICS 

Here we prove the robustness of the quantized velocity 
plateau dynamics against thermal fluctuations, changes 
in vertical load and the effect of quenched disorder in 
the confining substrates. In our simulations, we take 
into account complete layers, realizing an essentially crys- 
talline lubricant configuration at the given temperature, 
generally kept below its melting temperature. To in- 
vestigate the dragging of kinks and the ensuing exact 
velocity-quantization phenomenon, we evaluate the ra- 
tio w = v cmx /v cx t of the time-averaged lubricant CM 
sliding velocity to the externally imposed sliding speed 
v cx t stays pinned to an exact geometrically determined 
plateau value, while v ox t itself or temperature T or even 
relative substrate corrugation amplitudes are made vary 
over wide ranges. In detail, the plateau velocity ratio 
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is a function uniquely determined by the excess linear 
density of lubricant atoms with respect to that of the 
bottom substrate, thus of the length ratio Ab, Eq. ([3]). 
Despite not affecting the w p i a t value, the top length ratio 
At also plays an important role, since it sets the kink cov- 
erage 6 = Akink/^t = (l — A^ 1 ) At: assuming that the 
ID mapping to the FK model sketched in Ref. makes 
sense also in the present richer 2D geometry, the cover- 
age ratio should affect the pinning strength of kinks to 
the top corrugation, thus the robustness of the velocity 
plateau, as we will show in Sec. IIVI below. 
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FIG. 2: (color online). The time-averaged velocity ratio 
w — v clnx /vcxt of a single lubricant layer as a function of 
the adiabatically increased top-slider velocity Vext for differ- 
ent temperatures of the Langevin thermostat (panel (a)) and 
for distinct degrees of quenched disorder in the bottom sub- 
strate (panel (b)); atomic random displacements are taken 
in a uniform distribution in the interval [— E/2, S/2] horizon- 
tally and [—E/4, E/4] vertically away from the ideal positions 
of Eq. ((5j). All simulations are carried out with a model com- 
posed by 4, 29 and 25 atoms in the top lubricant and bot- 
tom layers respectively, with an applied load F = 25. The 
plateau velocity ratio (dot-dashed line) is w p i a t = 55 — 0.138, 
Eq. (H2J). 



As a convenient system for practical calculations we 
consider a bottom substrate made of 25 particles, and 29 
lubricant particles in each layer, i.e. Ab = 29/25 = 1.16, 
which produces 4 kinks every 29 lubricant particles in 
each layer. This value of Ab is not to be considered in 
any way special: we found perfect plateau sliding for 
many other values of Ab- Moreover we also investigated 
the plateau dynamics for an anti-kink configuration Ab = 
21/25 = 0.84. 

Figure[^a) reports the time-averaged horizontal veloc- 
ity Vcmx of the single-layer lubricant CM, as a function 
of the velocity w ox t of a fully commensurate top layer 
(0 = 1) for several temperatures set by the Langevin 
thermostat. The velocity ratio w = v cmx /v cxt is gener- 
ally a nontrivial function of w GX t , that for low temperature 



(here T = and T = 0.1) displays a wide flat plateau fol- 
lowed by a regime of continuous evolution. The plateau 
extends over a wide range of external driving velocities, 
up to a critical depinning velocity w C rit> whose precise 
value is obtained by ramping i> ext "adiabatically" . Be- 
yond t>crit, the lubricant leaves the plateau speed and, 
at finite temperature, moves up toward the symmetrical 
speed w = 0.5, as dictated by the thermostat dissipation. 
At zero temperature, the plateau state is abandoned with 
a discontinuous jump of w. The small-w oxt side of the 
plateau is difficult to address by numerical simulations, 
since the relative uncertainty in the determination of w 
increases due to w cxt -independent thermal fluctuations of 
(t) . By running longer averaging simulations at low 
Uexti we mitigate the error bars over w, and the data are 
consistent with a plateau dynamics extending all the way 
down to the static limit u cx t — > 0, like in the ID model. 5 
In particular, within the plateau, the T — calculation 
shows perfect matching of the mean lubricant velocity 
to the geometric ratio w p i a t of Eq. (fT2")) . By increasing 
the temperature, w deviates more and more from the 
perfect plateau value. At the highest temperature con- 
sidered, ksT = 0.5 e, near melting of the LJ crystal at 
fc B T m ~ 0.7er^ and even at the smaller k B T = 0.25 e 
no strict plateau is observed in the simulations. Note 
however that for w cxt < 0.5 the plateau attractor tends 
to affect the dynamics even in the presence of large ther- 
mal fluctuations. We have checked that finite-size scaling 
shows essentially no size effect on the plateau, and in par- 
ticular on its boundary edge w cr it- 

The dynamical pinning mechanism of the present 2D 
model is similar to that discussed for the simpler ID 
model. 5,5 The bottom layer produces a corrugated po- 
tential energy which, with its near-matching periodicity, 
is responsible for the creation of kinks, generally con- 
sisting of weak local compressions of the lubricant chain, 
with two lubricant atoms trapped in the attractive region 
between two bottom-substrate atoms. In the quantized- 
velocity state, these kinks become pinned to the upper 
slider (whose potential varies slowly on the lubricant in- 
terparticle scale) and are therefore dragged along. 

We probe the robustness of the quantized plateau state 
even in the presence of an "irregular" (i.e., quenched- 
disordered) bottom substrate, by running simulations 
with its atoms randomly displaced, horizontally and ver- 
tically, away from the regular positions, Eq. ([5]), by 
a significant fraction of the lattice equilibrium spacing 
ab- Calculations still display the perfect plateau ve- 
locity, which however terminates at a lower v cl -n. Fig- 
ure 03 shows the results of one specific realization of the 
quenched disorder of 3 different amplitudes S. By av- 
eraging over different realizations of disorder with the 
same amplitude £ = 0.2 ab, we obtain, for example, 
Wcrit = 1.47 ± 0.27, compared to the value v cr i t — 2.3 
obtained for the case a perfectly crystalline bottom wall 
at T = 0. 

Our overall evidence is that, in the 2D system as much 
as in the more idealized ID model, the phenomenon of 
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FIG. 3: (color online). The critical depinning velocity for the 
breakdown of the quantized plateau as a function of the ap- 
plied load per top-layer atom F. For the 1-layer curves, the 
same geometry as in Fig. [2j temperatures ksT = 0.01 (cir- 
cles) and 0.25 (diamonds) are considered. The 2-layer curve 
(squares) is computed with a doubled number of lubricant 
atoms (58 rather than 29) confined between the same sub- 
strates in the same z-range. 
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FIG. 4: (color online). Variation of the plateau boundary ve- 
locity Wcrit as a function of the kink coverage O -1 = iVt/iVkink 
for a lubricant monolayer and bilayer. Calculations show lo- 
cal maxima of v CI i t for commensurate values of O for both 
Mayer = 1 and 2; except at kink coverage 0=1 n cr it is gen- 
erally larger for Mayer = 1 than iVi ayor = 2. Simulations are 
carried out with F = 25, T = 0.01, and A b = §§. 



IV. SOLITON COVERAGE AND LUBRICANT 
SOFTNESS 



velocity quantization is stable, reproducible and ubiqui- 
tous, and not the result of some careful tuning of pa- 
rameters. In particular, Fig. [3] shows that the depinning 
velocity marking the end of the quantized plateau is a 
smooth function of the load F applied to the top layer 
over two orders of magnitude. Calculations suggest that 
a moderate load near unity is beneficial to the quantized 
sliding, at least at the low temperature ksT = 0.0 le. 
This is the result of a competition between the bene- 
ficial role of load that limits thermally induced "slips" 
of the kinks, and the detrimental effect of limiting the 
vertical lubricant movements at high loads. Indeed, at 
higher temperature, where thermal fluctuations are more 
active in destabilizing the quantized state (thus yielding 
a smaller u C rit)> calculations show a general slow increase 
of u cr it with F, favored by the suppression of thermally 
induced "slips" promoted by the stronger load-induced 
confinement. 



Figure [3] shows that in the large-load low-temperature 
regime, the critical velocity for Mayer = 2 is even larger 
than for iVi ayer = 1. In this plateau regime, we identify 
"horizontal" kinks in the lubricant layer adjacent to the 
bottom potential, while the upper lubricant layer shows 
weaker x-spacing modulations. 



The perfect matching of the number of kinks to the 
number of top-atoms 8 = Mmk/M = 1, as considered 
in all previous figures, is an especially favorable circum- 
stance for kink dragging, thus for the quantized-plateau 
phenomenon, but not one to be expected to occur easily 
in actual lubricated sliding. It is therefore important to 
examine situations where this kind of full commensura- 
tion is absent. 

In the following we study the behavior of the depinning 
point ii cr it as a function of the commensuration ratio 0. 
In particular we keep Ab fixed (i.e. the density of soli- 
tons), and vary the number of surface atoms in the top 
substrate. In Fig. [4] the depinning velocity v CI it, evalu- 
ated by adiabatically increasing of v ex t, is reported as a 
function of the commensuration ratio _1 = Nt/Nkiak- 
While the dependence on kink coverage is rather erratic, 
we can recognize that v cr n tends to peak at or near inte- 
ger values of O. In contrast, a fractional displaced from 
commensurate regimes produces a weakening or even the 
loss of the quantized plateau, represented by v cr it — 0. 
Figure [4] shows that this drop in u cr it is especially sharp 
for Mayor = 2, for which many incommensurate cover- 
ages show no quantized plateau, even at the smallest 
v ex t accessible within practical simulation times. We con- 
clude that the occurrence of quantized lubricant velocity 
plateaus, although comparably robust and widespread, 
may not be trivial to observe in practice: strongly incom- 
mensurate coverages and/or thick lubricant films could 
prevent its manifestation. 

The detailed dependence of the depinning transition 



FIG. 5: (color online). Variation of the plateau boundary ve- 
locity v cr it as a function of the kink coverage O -1 = A^/AWic 
for a lubricant mono-layer showing different interaction ener- 
gies with the substrates. Simulations are carried out with the 
same conditions as in Fig. [4] but with varied eu = £bi, (but 
unchanged interparticle interaction en = e). 
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FIG. 6: (color online). Snapshot of the dynamically pinned 
state of the 2D lubricant, showing the bond lengths as a func- 
tion of the bond horizontal position along the slider. Shorter 
bonds ~ <ib identify kink regions, while longer bonds w 1.16eib 
belong to in-register regions. Larger interaction with the sub- 
strates eti = 6bi > 1 favors the commensurate in-register re- 
gion. 



on coverage is affected by the relative "softness" of the 
lubricant. Figure O shows that when the interaction en- 
ergy eti — e bi of the lubricant with the substrates is in- 
creased to become larger than the intra-lubricant interac- 
tion en , as is the case for a soft lubricant wetting strongly 
the sliding surfaces, the quantized-plateau region extends 
generally to larger w oxt and to a wider range of coverages 
0. The reason for this is illustrated by the configuration 
snapshots of Fig. [5] a comparably stronger interaction 
with the substrates eti = £bl > 1 favors the commensu- 
rate in-register lubricant regions, shrinking the kink size. 



layer 

FIG. 7: (color online). Critical depinning velocity v CI i t as 
a function of the numbers jVi aycr of lubricant layers. All 
simulations are carried out in a condition that favors the 
quantized-velocity sliding state: the model is composed by 
4, 29 ■ Afi aycr and 25 atoms in the top lubricant and bottom 
layers respectively (thus = 1), with an applied load F — 25 
and T = 0.01. The data show an optimal dynamical pinning 
at ./Vi aycr = 3 and a tendency for v CI it to drop considerably 
as the lubricant thickness increases beyond that value. For 
A^iaycr > 7 no quantized plateau could be detected. 



Localized kinks are more suitable to being picked up and 
dragged by the corrugations of the top substrate. 

Figure [7] shows how the lubricant thickness affects the 
dynamically pinned state, i.e. i> cr it as a function of the 
number Mayer of lubricant layers in the fully commensu- 
rate = 1 condition. The data show quantized-velocity 
plateaus, displaying a maximum extension for Mayor = 3, 
followed by a decaying w cr i t as the lubricant thickness 
increases beyond that value. For Mayor > 7 layers we 
could detect no velocity plateau, at least within the v ex t 
range accessible to a practical numerical integration of 
the equations of motion. Vertical corrugations of the lu- 
bricant induced by the kinks originating and propagating 
from the bottom substrate have the effect of mediating 
the kink tendency to pin to the top-layer corrugations, 
therefore favoring the observed perfect velocity quanti- 
zation. For increasing Mayor, these z-displacements de- 
cay rapidly, and this mechanism becomes less effective to 
support the dynamically pinned state. 

Our calculations show that while the quantized plateau 
occurs when the solitons - misfit dislocations - are driven 
gently through the solid lubricant, the speed- induced 
transition to the u cx t > u C rit dynamically unpinned state 
can involve the melting of the lubricant, due to the large 
Joule heating induced by shearing. By analyzing the de- 
tail of the depinning transition for 0=1, we observe 
two different depinning behaviors to the left and to the 
right, respectively, of the maximum shown in Fig. [71 For 
-Mayer < 3, the depinning transition changes the solid 
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lubricant film into a disordered diffusive liquid, with a 
lower density, (thus characterized by an increase in rt7) . 
By contrast, a lubricant composed of four or more layers 
remains essentially solid through the depinning transi- 
tion. This is the case because the Joule heating pro- 
duced by shearing is removed with sufficient efficiency 
by the thermostat so that the effective lubricant tem- 
perature remains moderate. The lubricant remains solid 
in particular when u cr ;t is small, as is the case for thick 
lubricant layers, but also for 7^ 1, or if a larger damp- 
ing 77 ~ 1 is introduced in the Langevin dynamics. This 
indicates that the quantitative details of the findings of 
the present section depend on the assumed dissipation 
mechanism, and could therefore vary between one or an- 
other specific experimental situation. Nonetheless the 
qualitative result that under suitable conditions the de- 
pinning could be associated to lubricant melting while 
under other conditions the lubricant could remain solid 
through the transition is expected to be real, and should 
be verifiable in real boundary-friction experiments. 



V. FRICTION 

The (x-directed) kinetic friction force Fk that must 
be applied to the top slider to maintain its sliding mo- 
tion balances instantaneously the force that the lubricant 
exerts on the top substrate. Fk can be decomposed in 
several terms by analyzing the mechanical work done by 
the individual force contributions. The forces that the 
sliding top exerts on the lubricant [first term in Eq. ^] 
produce a total work VFt_ >i whose mean contribution to 
dissipation can be computed by evaluating the change in 
lubricant kinetic energy over a time interval At 



E k (t + At)-E k (t) 



(13) 



E 



Fj + /damp j + fj(t) ■ fj dt' 
, 477. 



pt+At 

Fj- fj dt' - —E k At 



VN\ / v cm ■r t dt' + J2 / fj ■ fj dt'. 
Jt 4 Jt 



rt+At 



The first term in Eq. (| 14[) represents precisely the contri- 
bution of the LJ top-lubricant interactions to Wt— a, plus 
bounded terms which do fluctuate, contributing a limited 
[order 0(At) a ] potential-energy change in the limit of a 
long time integration. The second term is proportional 
to the kinetic energy averaged over the At time inter- 
val. The last term in (|14[) describes the correlation of the 
random forces and the lubricant velocities. By careful in- 
tegration over the trajectory produced by the Langevin 
dynamics, averaging over the distribution of the random 
variables fj (indicated by angular brackets), and recall- 



ing Eq. ([9]), this term is evaluated as 

rt+At 

fj(t')-fj(t')dt') = 



1 



t+At 



dt' / dt"(fj(t')-fj(t")) 



(15) 



(16) 



t+At 



= - / dt' 



8ijk B T 



t 

t+At 1 

dt'~ 



dt" ArjkBTSit' 
4r)k B T 



t") x 2 (17) 



At. 



(18) 



The factor 2 in line (|17p comes from summing x and z 
components, and the factor \ in the integral of line (fl~8"|) 
is the result of integrating the Dirac delta at the inte- 
gration boundary. We observe that, like the conservative 
potential energy terms, the lubricant kinetic energy de- 
viation E k (t + At) — E k (t) is of order O(At) over a long 
time integration. By combining Eq. (|14[) with the result 
of Eq. (fTSf . in this large- At limit the average dissipated 
work becomes 



El 



Ni — 
N\k B T - —mv c 



■ n 



—to [r t - v c 



0-n 



At 



(19) 
O(Af) . (20) 



The fourth and last term in square brackets is added to 
account for the work done by the top substrate directly 
against the dissipation forces of Eq. (fTTj) . We omit the ex- 
plicit indication of the averaging over the Langevin noise, 
which becomes irrelevant in the limit of a long averaging 
time At. 

This total dissipated energy is conveniently written in 
a form where we measure the particles velocities relative 
to the instantaneous center-mass velocity w cm . Indeed, 
by using ft — (ft — v cm ) + v cm and omitting corrections 
of order 0(At°), we can rewrite the dissipated energy as 



WUi = — {Ek, 

TO 



Nik B T) At 



(21) 



+ (f t 



At. 



where terms linear in (fi — v cm ) vanish due to the defi- 
nition of Vcm, and Ek cm is the kinetic energy in the in- 
stantaneous CM frame: 



TO X ^ 
~9~ 2-^, 



(22) 



The work that the top layer does on the lubricant is 
supplied by the kinetic friction force Fk [the top-substrate 
vertical motion only yields a work of order O(At )], and 
equals Fk v cxt At. This relation provides an instructive 
decomposition of the average kinetic friction force: 

Fk = —^W t ^ x (23) 



Vext At 

477 
+ rjN 1 

Vcxt 



E k . 



[ft 
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FIG. 8: (color online). Results of the same model as in Fig. [2] 
(F = 25, T = 0.01). As a function of the top-layer velocity 
v cxt adiabatically increased (circles) or decreased (squares) 
the three panels report: (a) the average velocity ratio w — 
"cm^/woxt, compared to the plateau value w p i a t = gg — 0.138, 
Eq. (I12|l ; (b) the average friction force experienced by the top 
layer; (c) the effective lubricant temperature, computed using 
the average kinetic energy in the frame of reference of the 
instantaneous lubricant center of mass, Eq. (|24l) . 



Equation (|23|) expresses the kinetic friction force Fk as 
the sum of two contributions. The first and most im- 
portant term is the contribution of the fluctuations of 
the lubricant-particle velocities in the instantaneous CM 
frame, reduced by the thermostat-set value. Even at 
T = 0, this term would vanish in the event that the 
lubricant moved as a rigid whole. The term in the final 
line can be interpreted as the friction force that dissipa- 
tive phenomena with the bottom and the top respectively 
would produce, again even if the lubricant was rigid. This 
trivial term related to the overall CM motion yields a 
minimum total friction force whenever v cm = hf\, which 
corresponds to v cmx = h v cxt for the main horizontal 
component. In fact, this trivial contribution originates 
the tendency of w = v cmx /v cxt to abandon the plateau 
and reach 0.5, as in Figs.HJ and [HI [SI in the large sliding 
speed limit where this trivial term dominates. Equa- 
tion ([23]) can be reformulated in terms of an effective 
temperature of the steady state: 



B-icff 



1 

~Ni 



(24) 



^F\- 



The formulation (|24|) , clarifies the contribution of friction 
to Joule heating, and shows explicitly the irrelevant role 
of the trivial dissipation terms in this respect. 

Our MD simulations allow us to calculate instanta- 
neously the (x-directed) kinetic friction by considering 



FIG. 9: (color online). Results of a model composed by 4, 
21 and 25 atoms in the top, lubricant and bottom chains, 
Ab = 21/25 = 0.84, which according to Eq. ((12}, produces 
perfectly quantized dynamics at a negative u> p i a t = — 4r — 
—0.190, a dot-dashed line in panel a: this backward lubricant 
motion is caused by forward-dragged anti-kinks. The other 
simulation parameters are F — 25, T = 0.01. As a function 
of the top-layer velocity v cxt adiabatically increased (circles) 
or decreased (squares) the three panels report: (a) the av- 
erage velocity ratio w = v cnl x /v ex t; (b) the average friction 
force experienced by the top layer; (c) the effective lubricant 
temperature, Eq. 



the longitudinal force component that the lubricant ex- 
erts on the top driven substrate. By time averaging, the 
Fk value obtained coincides precisely with that evaluated 
through Eq. |23|) . 



A direct numerical evaluation of the friction force is 
shown in panels (b) of Figs. [8] and [9] for the kink and 
the anti-kink geometries, respectively, depicting typical 
scans where v ox t is cycled adiabatically up and down in 
small steps. Like in the T = simulations of the ID 
sliding model,— and as is well known for the static 
pinning/dcpinning of the standard FK model;2Ir— the 
present 2D model displays a clear hysteretic behavior 
of the velocity-plateau state (panels (a)), provided that 
thermal fluctuations are much smaller than the energy 
barrier hindering the full exploration of phase space for 
any affordable simulation time. A clear correlation of Fk 
with the effective lubricant temperature feeTcS oc Ekcm 
(panels Fig. [8h and Efc) is evident, but the trivial con- 
tribution is also seen to play a quantitative role. The 
bistability region with a clear hysteretic loop allows us 
to gauge the tribological effect of the quantized plateau 
state. In the kink geometry, Fig. [8] makes it apparent 
that friction is significantly less in the pinned quantized 
state than in the depinned state, all other parameters 
being equal. Figure [9] shows a similar tribological behav- 
ior, but with the occurrence of an additional interme- 
diate approximately quantized velocity plateau (the ID 
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0.25 




0.15 



FIG. 10: (color online). The average velocity ratio w — 
«cmi/»cxt as a function of the adiabatically increased (circles) 
or decreased (squares) top-layer velocity tw. The model pa- 
rameters are the same as in Fig. [8] but with no thermal effects 
(T = 0). 



model displayed similar intermediate plateaus 5 ). The dif- 
ference in friction between the primary and the secondary 
plateaus is a small one, with slightly less friction in the 
latter state. Comparison of the curves in panels and 
demonstrates that this effect is entirely due to the triv- 
ial friction term. This term is very large in the primary 
plateau because of the great difference of w from the 
optimal symmetrical drift w = 0.5; in passing to the sec- 
ondary plateau, the trivial friction term decreases quite 
substantially, but the kinetic fluctuation contribution to 
friction increases by a similar amount. 



VI. INTERMITTENT DYNAMICS AT THE 
PLATEAU STATE 



So far, we assumed some fixed externally imposed 
speed for the top slider. In a standard tribological sim- 
ulation, the slider is instead generally pulled through a 
spring of given stiffness whose other end is moved at con- 
stant velocity^ The spring can be viewed as a way to 
mimic not only the experimental driving device, but also 
the elasticity of the sliding substrates. From an experi- 
mental point of view, frictional forces may display a typi- 
cal sawtooth dependence at sufficiently low driving exter- 
nal velocities, the hallmark of the intermittent stick-slip 
dynamics. The details of this low-driving behavior of 
course depend sensitively on the mechanical properties 
of the device that applies the stress. 

Beside constant speed, we also explored the model with 
the spring pulling method. We obtained the quantized 
sliding state, also with the spring, and the simulations 
showed differences with the constant-speed method only 
in the fine detail. For this reason we omit here to display 
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FIG. 11: (color online). The time evolution of the average 
velocity ratio w and the corresponding kinetic friction Fk for 
the four dynamical states (a), (b), (c), and (d) marked in 
Fig. HQ] The first three panels, referring to quantized sliding 
states, display a typical intermittent stick-slip dynamics with 
small amplitude fluctuations; the last panel, not correspond- 
ing to a plateau state, exhibits chaotic large jumps in both w 
and Fk- the dot-dashed lines highlight the quantized-plateau 
value Wpiat- 



and analyze those results in detail. 

In fact, even within the rigid-drive model, for suffi- 
ciently low sliding speed u ex t and for extremely small 
thermal fluctuations (i.e., T ~ 0), we did observe a stick- 
slip intermittent dynamical regime in the total force ap- 
plied to the sliding top substrate in order to keep its slid- 
ing speed constant. Figure ITD1 depicts a typical adiabatic 
up-down scan of f cx t, similar to Fig. [8^ but at T = 0. 
To get insight into the different dynamical regimes out- 
side and inside the bistability region, Fig. II II analyzes in 
detail the four points marked in Fig. 1101 three of which 
belong to the plateau state and one outside of it. Each 
of the four panels of Fig. [TT] compares the time evolu- 
tion of the lubricant CM velocity to the corresponding 
instantaneous kinetic friction force Fk- 

For the first two plateau points (a) and (b) at low v cx t , 
the typical sawtooth time dependency of Fk is clearly 
visible. In these regimes of motion, the rescaled lubri- 
cant CM velocity w = v cmx /v ex t performs tiny periodic 
oscillations around the quantized plateau value and Fk 
exhibits an intermittent, almost regular, stick-slip pat- 
tern. This oscillatory phenomenology is understood as 
the kink array being dragged at full velocity w ox t across 
the associated periodic Peierls-Nabarro potential. 24 In- 
deed, the Fk oscillation frequency matches the washboard 
frequency of the kinks Vcxt/ab multiplied by their num- 
ber iVkink = N t = 4 in these = 1 simulations, due 
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to the four equally-spaced kinks crossing sequentially the 
corresponding Peierls-Nabarro barriers. A much less reg- 
ular pattern would arise in a 8 ^ 1 geometry. Due to 
the rather extended nature of solitons and smaller am- 
plitude of the Peierls-Nabarro barrier with respect to 
the full atomic corrugation, the quantized plateau stick- 
slip regime is associated to much lower dissipation than 
regular stick-slip dynamics as seen in dry friction AFM 
experiments! 25 ' 26 By further increasing the external driv- 
ing speed to v cxt = 1.8, panel (c), the kinetic friction 
becomes a little more "erratic", with an extra modula- 
tion of the sawtooth shape, resembling a sort of inverted 
stick-slip. Panel (d) displays the dynamics for the same 
v cx t = 1-8, but obtained while cycling i> oxt down: this 
dcpinned state shows large fluctuations both in w and 
.Ffc, giving rise to significant tribological dissipation, as 
expressed by the first term of Eq. ([23)) . 

VII. DISCUSSION AND CONCLUSIONS 

The present simulation work is meant to character- 
ize tribologically the quantized sliding state discovered 
in a ID sliding model 5 - and later observed in a substan- 
tially less idealized 2D model." The perfectly quantized 
plateau is demonstrated here to extend over broad pa- 
rameter ranges of the model, being robust against the 
effects of thermal fluctuations, quenched disorder in the 
confining substrates, the presence of confined multiple 
(up to 6) lubricant layers, and over a wide interval of 
loading forces. When temperature becomes comparable 
to the lubricant melting point, the plateau state tends 
to deteriorate and eventually disappears, but the geo- 
metrically determined velocity ratio u> p i a t, acting as an 
attractor, leaves anyway a trace in the ensuing "noisy" 
dynamics. 

The velocity plateau, as a function of w ox t, ends at a 
critical velocity w cr it, and for u ext > w cr ; t the lubricant 
tends to accelerate toward a speed 0.5v cxt - This result, 
dictated by the symmetric choice of the Langevin thermo- 
stat dissipation, may change in real systems where heat 
dissipation occurs through generally asymmetric sliders. 

Our calculations show that while the quantized plateau 
occurs when the solitons are driven gently through the 
solid lubricant, the speed- induced transition at v ox t > 
fcrit can involve the melting of the lubricant, due to 
the large Joule heating induced by large velocity shear- 
ing. The depinning value t> cr it is linked to the rate of corn- 
mensuration O of kinks to the upper slider periodicity: 



a particularly robust plateau signaled by a local maxi- 
mum of u cr it is located at well-commensurate O values, 
especially 6 = 1. Besides, the lubricant softness, set- 
ting the width of the propagating solitonic structures, 
is found to play a major role in promoting in-registry 
contact regions beneficial to this quantized sliding. Our 
typical dynamical depinning speed v cr it , is of the order of 
a few model units (corresponding to sliding speed values 
ranging from tens to thousands m/s for realistic choices 
of model parameters) , is very large compared to typical 
sliding velocities investigated in experiments. This sug- 
gests that in practice sliding at a dynamically quantized 
velocity is likely to be extremely robust. In real exper- 
imental systems, deterioration of the quantized sliding 
state will most probably associated to mechanisms such 
as disorder, boundary effects, or unfavorable lubricant- 
substrate incommensurate geometries, rather than to ex- 
cessive driving speed. 

By cycling v ext in underdamped regime, the layer slid- 
ing velocity exhibits a hysteretic loop around t> cr it , like in 
the ID models The bistability region allows us to gauge 
the tribological effect of the quantized plateau state. So, 
in the framework of Langevin dynamics, by evaluating 
the force instantaneously exerted on the top plate, we 
find that this quantized sliding represents a dynamical 
"pinned" state, characterized by significantly low values 
of the kinetic friction. A characteristic backward lubri- 
cant motion produced by the presence of "anti-kinks" , 
has also been observed in this context. This peculiar 
backward motion is likely to represent the most curious 
evidence of the quantized plateau state when - as we 
hope - it will be investigated experimentally. 

On the theoretical side, the role of substrate deforma- 
bility and defects in the lubricant structure, together 
with a realistic three dimensional description, with force 
fields representative of a concrete lubricated configura- 
tions, will certainly require further investigations. 
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